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With the advent of the Mariner ’71 Mission, NASA has been sending spacecraft 
to orbit various distant bodies within the solar system. At present , there is still no 
adequate theory describing the inherent state estimation accuracy, based on two- 
way, coherent range-rate data. It is the purpose of this article to lay the groundwork 
for a general elliptic theory, and in addition to provide an analytic solution for the 
special case of circular orbits. 


I. Introduction 

When one begins to analyze a possible future mission 
employing an orbiting spacecraft, it is desirable to know 
the characteristics of the position and velocity errors of 
that spacecraft, based on a single revolution of range-rate 
data. To know these errors as a function of orbit shape, 
size, and orientation with respect to the Earth is the cen- 
tral issue. If range were the data type of interest, the prob- 
lem would be vastly simpler, since velocities would not be 
involved. Unfortunately for the analyst, however, range- 
rate data is central and must be treated. This particular 
problem has proved to be a most formidable analytic en- 
deavor. 


II. Development of the Observable 

One of the most important aspects of this orbital es- 
timation problem is to pose it properly. In this regard, the 
“binary star problem” of classical astronomy can provide 
insight. In attempting to obtain the orbits of a distant bi- 
nary star system through the use of spectroscopic Doppler 


shifts, astronomers discovered that they could determine 
all the elements of the motion except one. The one ele- 
ment, which turned out to be completely indeterminate, is 
the longitude of the ascending node relative to the plane- 
of-the-sky, with the plane-of-the-sky being that plane nor- 
mal to the line-of-sight from the observer to the binary 
system. 

The reason for this indeterminacy is that the binary 
star system is effectively “at infinity,” meaning that no 
significant parallactic motion exists between the observer 
and the observed. 

In this problem, the body about which the spacecraft 
is in orbit is not “at infinity,” but is relatively near. This 
implies that there will be some parallactic motion between 
the observer and the observed, and hence we should be 
able to obtain a “reasonable” solution for the longitude of 
the ascending node relative to the plane-of-the-sky. This 
parameter will typically be the principal contributor to the 
spacecraft state error. As such, classical elements relative 
to the plane-of-the-sky appear to be a sensible “eigen” 
system in which to solve the problem. 
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The geometry of a spacecraft in orbit about a distant body is depicted in Fig. 1. Let: 






( 1 ) 

( 2 ) 


where x p , y p , and z p are the plane-of-the-sky Cartesian position coordinates of the spacecraft and x p , y p , and z p are the 
plane-of-the-sky Cartesian velocity coordinates of the spacecraft. 

The orbital coordinates are then defined as: 



with 
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where E is the eccentric anomaly, a is the orbit semi-major axis (km), e is the orbit eccentricity (dimensionless), and 
n is the orbit mean motion (s -1 ). 

The orbital coordinates can be related to the Cartesian, plane-of-the-sky coordinates as follows: 


r p = R T (Sl,L>,i)r 0 


( 12 ) 


where f? T (fi,w,i) is the transpose of the Euler rotation matrix, defined as 
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sin i sin w \ 


R{Q,,<jj y i) 
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\ sin 12 sin i — cos 12 sin i cos i ) 


(13) 


and 12 is the longitude of ascending node relative to the plane-of-the-sky, w is the argument of periapsis relative to the 
plane-of-the-sky, and i is the inclination relative to the plane-of-the-sky. 

To have any hope of analytically solving this problem, all possible simplifications must be made without, of course, 
loss of generality. One such simplification will be to define the x-axis of the plane-of-the-sky system to be the line-of- 
nodes. This implies that 


£2 = 0 


(14) 


Thus R t ( 12,0), i ) becomes 




/ cos u> — sin u? 

sin cj cos i cos u> cos i 
\ sin LJ sin i cos c j sin i 


0 \ 
— sin i 

cos i / 


If it is now assumed that the orbital elements are constant in time, 


(15) 


r p = R T (u), i)r 0 (16) 

From Fig. 1, it is seen that 

r d = r p — r T (17) 

( Vjyt COS 9 \ 

Vjsft sin 6 (18) 

VRt + pet / 

where: 

Vs is the component of the velocity of the target central body relative to the Earth that lies in the plane-of-the-sky; 
this shall hereinafter be referred to as the “cross-velocity” 

Vr is the component of the velocity of the target central body relative to the Earth that lies along the line-of-sight 

Pet is the distance along the line-of-sight between the center of the Earth and the center of the target central body 

at the initial epoch of estimation, i.e., at t = 0 

9 is the angle between the cross-velocity and the line-of-nodes 
One can now form the range-rate observable, p , as 
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where 


p = M 


( 20 ) 


Therefore, 


where 


It can then be shown that 


pp = r p * r p - r p fr - ry • r p + r T • r T 


ir = 


( Vjy cos 0 
V/v sin 6 
V R 


t t • r T = 4- 4- VrPet 

r p ip — a 2 ne sin E 

r p • tt = Vat (^ ori 4- ^ 2 ) 4- Vr(£ sinw 4- 77 cos w) sin i 

r T -r P = Vtft((ai 4- »)of 2 ) 4- (Vh< 4- />£t)(£ sinu> 4- r)cosu>) sin i 

where 

ai = cos 0 cos w 4- sin 0 sinw cos i 
(*2 = — cos 0 sin u> 4- sin 0 cos u cos i 

Next one has 

P = (r d ■ r d y = [r p • r p + r T • r T - 2r p • r T ] 3 

This can be shown to be 

p — |a 2 ^1 - e cos E^j + p% T + 2 petVrI + + V^jt 2 

— 2V/v<(^ai + r?a 2 ) - 2 (Vr* + Pet)(Z sinw + » 7 Cosu>)smi 


( 21 ) 
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(25) 

(26) 

(27) 
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(29) 
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For large pet this implies 


1 
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1 
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2 sin i f 
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This then becomes 


1 1 

- « 

P Pet 


i _ Y*L + sint 

PET PET 


sin u + r) cos u) j* j 


Combining Eqs. (21), (23)— (26) , and (32) yields 


p « — — IpExViR — Pet\£ sinu> 4- ») cost*;} sin i + Vj$t + o^ne sin E 
Pet L l J 

~Vjv|(^ + £t)a i -h (77 + f}t)ot2^ ~ sinu; + fjcosu; j|£ sin a; + tjcosu; j sin 2 i 


(32) 
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From the above, one can observe that 


|^sinu; + r/cosc*;} = a|(cosi£ — e)sinu; + /?sin E cosu;| 

sinu; 4* 77 cos u>j = — e~s~£?) {*- sin i? sinu; 4 (3 cos E cos a; | 


(34) 

(35) 


Combining Eqs. (33)-(35) produces 


( an sin ^ \ f . _ . „ _ 1 

Vr — - — { — sin E sin u 4 p cos E cos u > 

\ 1 — e cos E J 1 J 

H — { V£t — ( SSL* ) {sin E sinu; — /?cos£ , cosu;}{(cos E — e)sinw 4/? sin E cosu>} 

PET { \l - ezosE J l J l J 

+ a 2 ne sin i? - Vv{(f + £<) a i + (*1 + jj 


(36) 


Now one proceeds to expand the last term of Eq. (36). To do this, Kepler’s Equation is needed: 


E — e sin E = n(t — Tp) 


(37) 


where Tp is the time of periapsis passage. 
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Without loss of generality, one may assume that Tp = 0. Therefore 


nt = E — e sin E 


( 38 ) 


Combining Eqs. (5), (6), (8), (9), and (38) produces 


{(£ + fO«i +(*? + ^)a 2 } 

r / „ (i? — esin#) sin ./ . _ (£ — e sin E) cos E \ x 

= a\a l lcosE-e-± } — +a 2 [sinE+i } — \ 

l\ 1 — e cos E / V 1 — e cos E ) > 


(39) 


The range-rate approximation then becomes 


p & \ Vr — ( qns ? n * j f gin ^sinw + (3 cos E cos w } -+ — { V^t + a 2 ne sin E 

[ \ 1 — e cos E J l J pet L 

— ns * n |sin j£sinu> — (3 cos E cosu; j|(cos E — e) sinu; 4- /?sin E cosu;| 


- Vjva{«i(' 


cos E — e — 


(E — esin 2?) si nE 
1 — e cos E 


) + c* 2 /? (si 


sin E + 


(£ — esin £■) cos# 
1 — e cos E 


)}}] 


(40) 


If one wished to compute an approximation to the actual received range-rate signal, ignoring the motion of the 
receiving station relative to the center of the Earth, Eq. (40) would be the expression to use. In this case, however, the 
range-rate approximation shall be used for the development of partial derivatives only, and the only partial derivatives of 
interest are with respect to the spacecraft state (classical orbital elements with respect to the plane-of-the-sky). Hence, 
the elements of interest are a, e, Tp, i, and the final element is 0. Note that 0 by our definition is equivalent to the 
longitude of the ascending node, Cl. 

Thus those terms in the range-rate approximation which do not contribute to the appropriate partial derivatives 
will hereinafter be ignored. Additionally, those terms which contribute negligible partial derivatives will also be ignored. 
As a result, the range-rate approximation which shall be used for the computation of the partial derivatives is 


/ ansint \ f . n . ^ ^ 1 

P « “ — s sin E sin u — p cos E cos u > 

\ 1 — e cos E J ^ > 


-(£){*■( 


cos E — e — 


(E — esin £')sin 


1 — e cos E 


+ot 2 0 (si 


sin E + 


(E — e sin E) cos E 
1 — e cos E 


)}] 


(41) 


The partial derivatives with respect to a, e, Tp, u, and i shall be computed from the first term in the range-rate 
approximation (since these dominate partials for the same parameters taken from the second term), whereas the partial 
derivative with respect to 6 will be computed from the second term (this is the only term with 6 variation). 
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III. Obtaining the Partial Derivatives 

Rewriting Eq. (41) produces 


r _i 1 . 
a 3/x 2 sin 

1 — e cos E 


7 j |sin E sin u — 0 cos £* cos 


_ f 

VPET/ \ 


(£ — e sin £*) sin £\ . / . _ (E — e sin E) cos E 

<*1 ( cos E — e — - ( sin E + - } — 

1 1 - e cos E ) \ 1 - e cos E 


)}1 


(42) 


To produce the partial of the observables with respect to the classical state of the spacecraft, the partials of the eccentric 
anomaly, E } with respect to the classical state are needed. They are 



3 (E — esin E) 
2a(l — e cos E) 

sin E 

(1 — ecos E) 
n 

(1 — ecos E) 


= 0 
= 0 
= 0 


Taking partial derivatives of the observable produces 




{— sin E smu + /?cos jEcosu;} 
(1 — e cos E) 


3 (E — esin E) {(cos E - e) sin a; + /? sin E cos uj] 
(1 — ecos E) 3 



an sin i 


— cos E {(cos E — e)cosu — /?sin 2£sinu;} 
(3(1 - ecos E) 2 


+ 


sin E {(cos E — e)sina; -f /? sin E cos u > } 
(1 — e cos E) 3 


(43) 
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(48) 
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(51) 
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. [ sin E sin u — (3 cos E cos u> 


(1 — e cos £) 


- (£) M 
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(£ — e sin E) cos E 
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(52) 

(53) 

(54) 


where 

(55) 

(56) 


c *3 = — sin 0 cos cj + cos 0 sin cj cos i 
a 4 = sin 0 sin u; + cos 0 cos u; cos i 


Again, if one is to get through this problem, one must simplify as much as possible or risk getting lost in the 
algebra. A careful examination of the partials detailed in Eqs. (49)-(54) reveals that the first term of Eq. (49), other 
than a constant of proportionality, is the same as Eq. (53). Therefore it may be eliminated from Eq. (49) by a change 
of variable. Let 


which implies that 


From Eq. (57), one has 



a + ( 2 a cot i)i 


(g ) = (£) it) + (I) (£) 



Therefore, 


dp\ / — 3nsini\ (£ - e sin E) { (cos E_ - e) sin u + /? sin E cos w) 
aa 7 J - V 2 ) (1-ecos E) 3 


(57) 


(58) 


(59) 

(60) 


(61) 


To obtain the correct partials in Eq. (58), one must carefully discern which of the variables in Eq. (57) are differentiable 
parameters and which are not; it is then clear that Eq. (61) represents a considerable simplification of the problem. 
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IV. Accumulating the Information Array 

The normal, Gaussian way by which to obtain the covariance of errors in a set of estimated parameters, T, is stated 
as follows: 


r = [f ~ 1 +a t wa^ 1 


(62) 


where f is the a priori error in the covariance, A is the vector partial derivative of the observable (at the time of 
observation) with respect to the parameters being estimated, and W is the weighting matrix. If it is assumed that the 
error in the measurement of the observable at time t{ is uncorrelated with the error in the measurement of the observable 
at time tj , for all i and j, then W becomes a diagonal matrix of dimension equal to the number of observables. From 
minimum- variance concepts, the diagonal elements of the weighting matrix may be shown to be just the inverse of the 
variance of error in the observable (it is assumed the error process of the observable is stationary and hence that the 
variance of error in the observable is constant and unchanging for every data point). 

If one assumes the existence of no a priori information, then the error covariance sought, assuming discrete obser- 
vations, may be formed as the inverse of the information array: 


r = 



(63) 


where X is the vector of estimated parameters (in this case this is the spacecraft state in classical elements relative to 
the plane-of-the-sky), which shall be ordered in the following way: 


— (pi) P2, P3y P4y Pb, = ( 9 , a 7 , i, e, 7>, 


and N is the number of data points (observations). 

In this case, to simplify the accumulation, continuous observations are assumed, and hence the sums involved will 
become integrals. For further simplicity, the accumulation shall be performed over a full orbit period in the time domain 
P (in the manner in which “real” data is taken, in equal increments of time), beginning at —P/2 and ending at +P/2. 
Since the observable and the partials are all formulated in terms of the eccentric anomaly, E, it is logical to perform the 
accumulation with this parameter as the independent variable. From Eq. (37), one has 


dt = v 1 ~ e cos ~l dE 

n 

Thus a general information array element, /*j, may be stated as 


(64) 


lij 


1 r ,et\ (M\ 

(a<r}) y_ T V dpi ) \ dpj ) n 


(65) 


where a is a scale factor chosen such that the information array, Uj , obtained continuously is the same as that formed 
discretely. 
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To ascertain the value of a, let X consist of only one parameter, with its attendant partial derivative being a 
constant, K. Then the discrete information is just 


/n = 



The continuous information, however, is 


/n = 



( 66 ) 


(67) 


where P is the orbit period (in seconds). 

Therefore, for the two informations to be equal, one must have 

*-(£) (68) 

This, however, implies that 

a = At (69) 


where At is the discrete sample rate (in seconds) . 

In the standard case of two-way, coherent range rate as the orbital data type, the discrete sample rate is 60 seconds. 
Thus the elements of the information array may now be stated using the continuous formulation as 


- _J_ r (°L) (°L) 

l '> ~ (alAt) J.Adpi) \dpjJ 


(l-ecosff) i£ , 


(70) 


For simplicity, scaled information array elements, I-j , shall be stated where 


* ■ (**) 


(71) 


Note that Eq. (71) implies that when the determinant of the full 6x6 information array is computed, to obtain a correct 
value, the result must be multiplied by 
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Similarly for the cofactors of dimension 5x5, the result must be multiplied by 


(^) s 

Since the information array is symmetric, only the upper diagonal portion of the array shall be stated. 

If one examines Eq. (42), it can be seen that the range-rate value obtained for a certain value of u> is the negative of 
that for u + 7r. The partials will also be of opposite signs. The information array elements will, however, be the same. 
Hence the information array must be periodic in u>, but go through two full periods as u goes from zero to 360 degrees. 

With the complicated nature of the partials, there is a need for a very specific table of integrals to perform the 
accumulation. Such a table is provided in the Appendices. With the use of the Appendices, the elements of the 
information array can be straightforwardly computed. 

Due to the complexity of the information array, any further simplification will be most welcome. We are only 
interested in obtaining the variance (standard deviation) of 0, which is the error in the line-of-nodes relative to the 
plane-of-the-sky. Any element in the inverted information array (the covariance), , can be stated as 


_ U-\) i+ iCo D tJ 

11 I D 


(72) 


where D is the determinant of the full 6x6 information array, i 7 , with the appropriate constant factors included, and 
CoDij is the cofactor of the ij th element of the information array, which is just the determinant of the remaining 5x5 
array obtained by removing row i and column j from the full 6x6 information array, with also the appropriate constant 
factors included. 

In this case, we are interested in Cn only. This is just 


Cn — Vg — 


Co Du 
D 


(73) 


If it is assumed that there are removable coefficients, a,-, in each of the ith rows and ith columns of and similar 
coefficients, 6;, coming from the 5x5 matrix whose determinant is the CoD n , then 


Cn — — 


CoD\A 

& y 


UU b l = 1 ( CoD' n \ 

n - =1 «. 2 d > ) 


where Coiy n and D? are the corresponding determinants with common factors removed. 


(74) 


One can now state the information array, i 7 , to compute U and CoD \\ . The information array, with coefficients 
removed, is shown below. 
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( 75 ) 


/;, = [{y + 1 - 2« + 2e ! } («I + o|) 

+ (?) {' #SW + ^ ln { 2^} - J (» - «’> - ? + + < - |<» - w 1 )} (“5 - “I) 


where 


(a^ + a 2 j = sin 2 0 + cos 2 6 cos 2 i 


(76) 


(a 2 — a|^ = 2 sin 2u> sin 6 cos 6 cos i - cos 2 u (sin 2 9 - cos 2 6 cos 2 


1 + /? ^+2e + ±(l- /?)(!- 3/?)} 


7 ' 2 = [{ W‘ S(t) - 1 ln ( 1 + t) - + 0) + (T^) } 

-^{(21n(l + 7) + 31n|^ r ^} + 2e + | 
x (sin 0 cos 2w — cos 0 cos i sin 2 w j | 

l' 13 = 0 sin 0 + j? 2 + ^ ^ ln 1 } } (sin 0 cos 2 w - cos 0 cos i sin 2w ) 





/u = I ( — I < ln 


1 + 0 


(77) 


H 


I' l5 = cos * cosi ' - ( 2 (V+eK ( l 1 + t /^ ) ( s * n ^ sin 2« + cos 0 cos j cos 2w) ] 

/' 6 = ^[ co8<,cos *‘ - (t 2 + ! + ? 111 { 2[n4 }) ( si 


sin 9 sin 2w -f cos 9 cos i cos 2u> J J 


(78) 

(79) 


e/9/ r l2(l + e)} +2 " W + W }“°‘' 

^ |(2 — e 2 ) ln 1 2 (\~+ ) } + - - 0) + 2e - e 4 } (sin0cos2u> - cos0 cos»'sin2u;jJ (80) 


(81) 

(82) 


I'-n = ~ ji !“(! + 7) + £(1 - e) - ^ + 20+ 3)} 


„ fe 2 „.. (2 — 5e 2 ) (l-e)(2 + e 2 ) 1 (1 + /J)(3 - 0)' 

- co82w ( 4^ 5( 7)+ 3 e 2 /3~ n(1 7) + +4p 



(83) 

(84) 

(85) 
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3c)/? 3 


7 ”~ (l2ev) [' + a + e)» 

r “ = [rh-iT^]" ,,2 “ 
/Sa = 1 "(lfl)“* 2 " 


sin 2 cj 


r, - (!-/?) rnn n r , 

34 “ e/?(l + /?) cos2w 


/35 = 0 


736 - (rrl) 


sin 2w 


V = * 

44 4/?3(l + /?) 


2(3+/?)+ o^ cos2w 


7 « - (4^) 


sin 2a; 




iu= 


1 

4/? 5 


2(2 + e 2 ) — e 2 cos 2u> 


r -1 

^56 - 


■^66 = 1 + 


(£0 


cos 2w 


( 86 ) 

(87) 

( 88 ) 

(89) 

(90) 

(91) 

(92) 

(93) 

(94) 

(95) 

(96) 

(97) 


where the removable coefficients, of course, turn out to be the coefficients of the individual partial derivatives, which are 


£ 

11 

0 

(98) 

V Pet / 

( — 3n sin *\ 

(99) 

— ( 2 ) 

03 = an cos i 

(100) 

a 4 = an sin i 

(101) 

a 5 = —an 2 sin i 

(102) 

ae = an sin i 

(103) 
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From these scale factors and Eq. (71), one can infer that the error in 9 may be written as 


<r e = <Tp 



Pet(p)H A<)’ 

/(e,w) 

a p 

. V«r(a)i0 0* . 

.g(0,i,e,u)' 


( 104 ) 


where /(e,o;) is the cofactor of I? n , and g(9, *,e,o;) is just the determinant, D* . 

V. Inverting the Information Array 

Since the interest is in a single element of the inversion, only those computations required in Eq. (74) shall be 
performed. Before plunging into the general elliptic case, the inversion for the special case of the circular orbit (i.e., 
e = 0) shall be demonstrated. 

A. Circular Case 

After much work, the elements of the information array (in the limit as e tends to zero) can be shown to be: 


1 12 


1 22 


*23 


/ 24 


*25 


*26 



(105) 

7 r 2 . 

= y sm« 

(106) 

= sin 9 -f ^ ^sin 9 cos 2a ; — cos 9 cos i sin 2a; j 

(107) 

= sin 9 — i ^sin 9 cos 2a; — cos 0 cos i sin 2a;^ 

(108) 

= cos 9 cos i — i ^sin 9 sin 2a> + cos 9 cos i cos 2a; ^ 

(109) 

= cos 9 cos i — i ^sin 9 sin 2a; + cos 0 cos i cos 2a; ^ 

(110) 

7T 2 1 

= — — - cos 2a; 
3 2 

(111) 

1 

= - cos 2a; 

(112) 

= 1 — - cos 2a> 
o 

(113) 

— — i sin 2a; 

(114) 

= — ^ sin 2a; 

(115) 

~ ^44 = ^55 = ^56 “ ^66 “ 1 

(116) 

= /85 = *36 - '45 = *46 = 0 

(117) 
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To perform the inversion in this special, circular case, it is necessary to remove the parameter u> from the estimation 
since the argument of periapsis is indeterminate. Additionally the parameter Tp is removed from the estimation since 
the partials are large and the same as those of w, only with the opposite sign. This rank reduction corresponds to 
assuming that the true anomaly of the spacecraft at epoch is perfectly known, which is realistic in comparison to the 
error in 9. With a large amount of algebraic manipulation, the four-parameter inversion yields an error in 9 of 


where 



Pet{p)*( A<)* 

/(«) 

(Tp 

V N (a)i(r 0* . 

g(9,i,ui). 


/(“ 0 = 


7T 2 , 1 „ 13 j 

1 + - cos 2w - — cos* 2w 

3 6 36 


(118) 


(119) 


,(* , i,„) = (2d£?di) [„< - ,» cos 2* - cos 5 V - £ + f cos 2* - f cos 5 2u.] 


( 120 ) 


It is now apparent from Eqs. (118) — (120) that the use of range-rate data for state estimation of a circular orbit about a 
distant planet can yield singular results. Singular state errors will result if 


9 = 90 deg, 270 deg, and/or 

(121) 

i = -90 deg, 90 deg, and/or 

(122) 

O 

II 

£ 

(123) 


This latter condition of singularity may occur when the spacecraft is orbiting a superior planet (at the beginning and end 
points of retrograde motion). These are important results, which are significantly at variance with orbit determination 
experience for highly eccentric trajectories, which appear well behaved and certainly not singular. Note that even with 
zero eccentricity, the results of Eqs. (119) and (120) imply a functional dependence upon u. This is an artifact of the 
definition of the coordinate system at time zero. If we assume w = 90 deg, then the maximum variation in the error 
in 9 for any other value of u can be shown to be 13.6 percent. Thus the error in 9 can be said to be only weakly 
dependent upon w. 


B. Elliptic Case 

Many thanks go to R. J. Muellerschoen of the Navigation Systems Section for his help in using MATHEMATICA 
to attempt the 6x6 inversion. Even with a large amount of available memory, the problem quickly consumed it all and 
more. Thus it appears unlikely that a “reasonable” functionality for / and g can be obtained in the general elliptic case. 
However, in having the information array, we have performed the accumulation in a general manner which will save 
much time. A large number of orbit-determination studies as a function of geometry relative to the plane-of-the-sky can 
be made quickly and cheaply. 
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VI. Results 

A. Circular Case 


Venus Orbiting Imaging Radar (VOIR), the progenitor of Magellan, was extensively studied many years ago. It 
was in these studies that the circular orbit singularities were first observed. The nominal orbit just happened to have 
conditions that were near-singular. As a result, orbit determination errors were obtained that were outrageously large 
(larger than the nominal orbit). This led to the use of several different computer programs to examine whether there were 
program errors or not. The two main programs used to investigate the state estimation were the Orbit Determination 
Program (ODP) and the Sequential Orbit Determination (SOD) Program. Both of these pieces of software yielded 
the same, poor orbit determination behavior and so confirmed the belief that we had inadvertently stumbled upon 
singularities in the orbit determination process for circular orbiters. The important parameters defining the nominal 
VOIR orbit were: 


Epoch 

Semi-major axis, a 
Eccentricity, e 
Inclination, i 
Argument of periapsis, u> 
Cross-velocity, Vn 
Target range, p E r 
Angle, 6 
Period, P 


March 17, 1982 
6552 km 
1.52 x 10' 4 

78.17 deg (relative to the plane-of-the-sky) 

170.7 deg (relative to the plane-of-the-sky) 

14.413 km/sec 

8.36765 x 10 7 km 

85 deg 

1.624 hr 


The computer solutions were assumed to have as a data base one full orbit of coherent, two-way range-rate data. A 
slight problem complicating exact comparison with the theory is that the computer programs employed assumed that 
when the spacecraft was not visible to the Earth tracking station, there would be no data available. This eliminated 
about 36 minutes of data (depending on the nominal inclination relative to the plane-of-the-sky) out of a total of 98 
minutes available. The theory, however, assumes that a full orbit of data is available no matter what. As a result the 
displayed errors are not in perfect agreement. The trends as functions of the angle 6 and the inclination i, however, are 
in excellent agreement and as seen in Figs. 2 and 3, clearly demonstrate the singular behavior of the circular orbit. The 
best place to compare theory and computer results is in Fig. 3 where the value of the inclination relative to the plane- 
of-the-sky is small. In this configuration, a full orbit of data would be available to both theory and computer results, 
yielding agreement to about 30 percent. As the inclination increases, the agreement between theory and computer 
results degrades as more data is occulted. The error parameter displayed in these figures is 


\<r x \ - r<Te - a<?9 


(124) 


B. Elliptic Case 

The orbit determination studies for the general elliptic case are yet to be published. 1 


1 For additional information on the initial studies which helped to delineate the circular orbit problems, see the following: J. Ellis 
and R. K. Russell, “Earth-Based Determination of a Near-Circular Orbit About a Distant Planet,” JPL Technical Memorandum 
391-406 (internal document), March 9, 1973; R. K. Russell, “Improved Orbit Determination of Low-Altitude, Nearly Circular 
Planetary Orbiters,” Engineering Memorandum 391-561 (internal document), April 25, 1974; R. K. Russell, “Analytic Orbit 
Determination, A General Methodology and a Successful Example,” JPL Engineering Memorandum 314-59 (internal document), 
August 9, 1976; and L. J. Wood “Orbit Determination Singularities in the Doppler Tracking of a Planetary Orbiter,” Journal of 
Guidance, Control, and Dynamics , volume 9, number 4, p. 485, July-August 1986. 
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VII. Conclusions 


It has been seen that circular orbits about distant planets may suffer singularities in over-all position error estima- 
tion. These singularities are due to orbit inclination, placement of the line-of-nodes, and insignificant cross-velocity at 
the start and end of retrograde motion when orbiting a superior planet. 

Even though these conclusions appear to yield poor state estimation, one should not be unduly alarmed inasmuch as 
the stated conditions for singularity are not maintained for extended periods during typical mission scenarios. However, 
mission analysts should be aware of these potential pitfalls and realize that spuriously large results for circular orbiters 
can be obtained and are not the result of incorrect assumptions or faulty software. 

The general elliptic problem appears so involved that analytic inversion at this time is just not feasible, and in any 
case the resulting expression for the position error would likely be so lengthy that any understanding would be lost in 
the maze. 
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z (TOWARD EARTH AT t = 0} 



Fig. 1. Spacecraft orbital geometry in the plane 
of the sky. 



Fig. 3. Norm of position error versus inclination, i, 
to the plane of the sky (with 0 ~ 0 deg). 




Appendix A 
Required Integrals I 

To perform the accumulation of the information array, a large number of difficult integrals of mixed trigonometric 
and algebraic functions and their powers are required. Contained herein are the necessary results. Note that due to the 
limits of integration, all integrals of odd functions will vanish. 

An important, general integral, from which some of the results contained herein are derived, is 




dE 




+ b cos E) n+1 ( a 2 _ 52)^ 


rPn 



(A-l) 


where P n are the Legendre polynomials. 


This integral is found in [1; p. 383, 3.661, 4]. 
For simplicity and convenience, we shall define: 


P — \J 1 — e 2 

The integrals necessary to perform the accumulation are: 


J (1 — e cos E)dE — 2 tt 
/> — e cos E) cos EdE = — 
J (1 — e cos E) cos 2 EdE — ? 

b°= 

L 

J cos 2 EdE = it 
J cos 3 EdE = 0 
J E sin EdE = 27 r 
J E sin E cos EdE = — j 


27T 

cos EdE = 0 


(A-2) 

(A-3) 


(A-4) 
(A-5) 
(A-6) 
(A-7) 
(A-8) 
(A-9) 
(A-10) 
(A-l 1) 
(A-12) 
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(A-13) 



dE / 2jt\ 

(1 — ecos E) \ /? ) 

cosEdE 

(1 — ecosE) \e/3j { P) \0 ) 
cos 2 EdE 


(1 — e cos E) 

cos 3 EdE 
(1 — ecos E) 

cos 4 EdE 
(1 — e cos E) 

E sin EdE 


=(*)< i -» , < 2+ »=(^) 

_/2»v r i+/? i 

(1 — ecos jE 7) \ e / n 1 2(1 + e) / 
EsinEcosEdE /2tt\ [" 1 ( 1+/? 11 

(1 — ecos £7) " V e 7 L * \2(l + e)/J 


T 2 (2 + /9) 


E'sinJE'cos 2 EdE 
(1 — e cos£) 


-&)[- 




i+/? 
2(1 + e) 


}]] 


(A-14) 
(A- 15) 
(A-16) 
(A-17) 
(A- 18) 
(A-19) 
(A-20) 


For a full development of the following integral, see Appendix B. 


r E 2 dE 

L t (1 — e cosl?) 



(A-21) 


where 


Continuing, 


5(7) = 


3 n 2 



E 2 cos 2 EdE f 2w\ [ ( 

(1— ,co»E) = (^jl 5(l,)+<; ( 2e -T)j 

dE _ /2i\ 

(1 — ecos E) 2 \/? 3 / 

cos EdE _ / 2ire \ 

(1 — ecosi?) 2 \ (3 3 ) 




(A-22) 

(A-23) 

(A-24) 

(A-25) 

(A-26) 
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£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 


= (£)*»+*« 


cos 3 

(1 — ecos l?) 2 

(!*£!£)» = (^ 5 ) (1 - « !(2 + v ^ - ^ = fe ) " ! < 2 + ^ - 2<J! - 

EsinEdE /2tt\ T 1 1 1 /2x\ , 

= (tJ b-i^J = UJ< 1 - T > 

1 + /? 


(1-ecos#) 2 

E sin E cos EdE 
(1 - ecos £) 2 

E sin E cos 2 EdE 


f?*\ [1 1 . f 1 + /? V 

~Uvb 1 + e + n \ 2(1 + e) J 


'2*\ n 

e 3 


1 


(l-ecosi?) 2 V e3 7L/? 1 + e 6 + { 2(1 + e ) } 


E 2 dE 


(1 — e cos E) 2 

dE 

(1 — e cos E) 3 

cos EdE 
(1 — e cos E) 3 

cos 2 EdE 


= (^) [s( 7 )- 4 / 31 n(l + 1 ) 




- ( Z7re \ 
~ U 5 ) 


(1-ecosE) 3 ~ U 5 J (3 ^ 


COS' 




(1 T e "g)» ■ ^ J (* ■ - + *» + 2 ^ - 2 « = (£) + + 2 ^ - 2 « 

■ (&) + * - 2 ^ ■ (&) + 2<J - *”> 

EsinEdE /ir\ T 1 1 

“ W [ P 3 ~ (1 + e) 2 . 

(1 + 2e)1 


(1 — e cos E ) 3 
i^siniJcos EdE 


(1 — ecos l?) 3 (e 2 ) [/? 3 ^ 2/? ^ + (l + e) 2 J 

Esin Ecos 2 EdE _ / * \ M_ _ 2) (3 + 4e) _ f 

(1 — e cos E) 3 ~ L 3 / [/? a( ** + (1 + e) 2 \ 

E 2 dE 


1 + P 


;}] 


(1 + e) 2 '“‘l2(l + e) 

(1 - e cos £)* = (?) [< 3 “ ^ - W1 " (1 + T) - 2e ' J<1 - T > 

E 2 cos EdE _ ^JLj | 3( 1 - - 4,3(3 - 2,3=) ln( 1 + - 2c/J( 1 - 

(l-«cos£)4 = (??) [(3-5^ + 2rtS(r)-4«3-4?)1n(l + 7)-2 e «l- T ) 


(1 — e cos i?) 3 

E 2 cos 2 EdE 


(A-27) 

(A-28) 

(A-29) 

(A-30) 

(A-31) 

(A-32) 

(A-33) 

(A-34) 

(A-35) 

(A-36) 

(A-37) 

(A-38) 

(A-39) 

(A-40) 

(A-41) 

(A-42) 

(A-43) 
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(A-44) 


h (5-4/? 2 ) 


r cos 3 EdE 


cos 4 EdE 


= (w>- 


(1 - /?) 2 ( 5 + 10/? - 10/? 3 - 4/? 4 + 2/? 5 ) 


= 7 2 (5 + 10/? - 10/? 3 - 4/? 4 + 2/? 5 ) 


f* EsmEdE / r ’r\[(3-^ 2 ) 2 

J_, (l-ecos^-UeJi /? 5 (1 + e) 3 

[* E sin E cos EdE _ (jr_\ [J_ _ 2 (. 

I (1-ecosE) 4 — V3e 2 / [/? 5 ^ ^ } + (: 


"Ge 2 )./? 5 ^ 4/?2 


(1 4- 3e) 
(1 + e) 3 


f E sin E cos EdE ( ^ \ ^ 7^2 , (*^4\ ^ /-i • 9 < o-2\ 

L (1 - ecosg) 4 °(3?) b (3 - 7 ' J +6 ' S )-(TTJ)5 (1+3e + 3e) 


r dE_ 

J_ w (1 — ecosE) 5 


(35 - 30/? 2 + 3Z? 4 ) 


(1-ecosE) 5 


(7 - 3/? 2 ) 


cos 2 EdE — ( JL 
(1— ecosE) 5 \4/?* 


(35 - 35/? 2 + 4Z? 4 ) 


cos 3 EdE ( 5ire 


/_" = (jf) < 35 - 40(,! + 8 ^> 


Esin EdE _ /w_\ r (5-3/? 2 ) _ 2 

(1 — ecosE) 5 \4e/ _ /? 7 (1 + e) 4 


(1 + e ) 4 


/ ' Esin Ecos EdE / ir \ T 1 /1C „.„2 , , 2 , . x 

(l-ec E )S = fe) IF 15 “ 2W + ^ + (TT^r + 4t) j 
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L 


* E sin E cos 2 EdE 


L 


(1 — ecos E)* 

vl2e 3 / 

E 2 dE 

fM[ 3 

(1 - ecos E)* 

Vi2/?v l 3 

— 2e/?( 1 - 

E 2 cos EdE 

f ^ [ 

(1 — e cos E ) 5 

\ 12e/? 9 / 1 
-2e/?(l - 

E 2 cos 2 EdE _ 

{ * ) 

(1 — ecos E) 5 

V 12e 2 /? 9 / 


(A-59) 


(A-60) 


(A-61) 


-2e/?(l - 7 )(57 + 12 0 - 66/? 2 - 8/J 3 + 12 /? 4 ) + 2e/? 2 ( 12 + 0 - 8/? 2 )] (A-62) 
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Appendix B 
Required Integrals II 

Out of all these difficult integrals, only one class is not expressible in closed form. These are integrals of the form 


E 2 cos m EdE 
_ x (1 — ecos E) n 


(B-l) 


We shall confine our interest to the lowest order of the above, namely where n = 1 and m = 0. The other integrals can 
be simply obtained from this one, which is 


-L-t 


E]dE_ 

ecosE 


(B-2) 


If the square function is expanded in a Fourier series over the interval — 7r to tt, one obtains 


E 2 = - 4 [cos E - cos 2E + ^ cos 3E j , —w < E < 7r 


(B-3) 


This now can be made to yield the integral we are seeking: 



E 2 dE 
1 — e cos E 


r — d JL — +4 f c ° sn£ ) dE 

3 y_ T l-ecos E ” 2 V 1 - e cos EJ 



cos nEdE 
1 — e cos E 


Thus, 



cos nEdE 
1 — e cos E 



cos nEdE 
1 — e cos E 


This latter integral can be found in [1; p. 366, 3.613] as 


cos nEdE 

-if 

i-/?V 

/ 0 1 — e cos E 

P V 

e ) 


(B-4) 

(B-5) 


(B-6) 


(B-7) 
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